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Abstract 

We discuss the predictions of perturbative QCD for angular flows 
of final state particles in two and three jet events including their cms 
energy and jet resolution (y cu t) dependence. The simple analytical 
formulae for gluon bremsstrahlung from primary partons, modified 
for gluon cascading, reproduce the main features of the experimental 
data well. For y cu t-selected events, the particle flow is derived from 
a superposition of colour dipoles in much the same way that photon 
radiation is derived from electric dipoles. 



1 Introduction 



We study the angular distribution of particles between the hadronic jets in 
the final state of a high energy collision. The angular flow of such particles, 
most of them with low momentum < 1 GeV, is expected to depend character- 
istically on the colour connections of the primary partons |H2|3j (for recent 
reviews, see, for example, Ref. 4J). An early example of such a phenomenon 
concerns the final state e + e~ — > qqg. The soft radiation is not distributed 
symmetrically between the three jets, but rather is depleted in the angular 
region between the q and q directions. This effect was first predicted within 
the string hadronization model [3] . Subsequently, it has been derived within 
perturbative QCD jB] where the particle flow is derived directly from the soft 
gluon bremsstrahlung, which is emitted coherently from all primary partons 
but with different strength from gluon and quark emitters according to the 
QCD colour factors. 

The first observation of the effect was due to the JADE collaboration [7j 
with many of the subsequent details having been studied by various experi- 
mental groups, for example |8p9pi0j . This "string/drag effect" is by now well 
established and reproduced by the popular Monte Carlo models. The ana- 
lytic calculations have been verified mainly for angles midway between the 
jets. Until now, the full angular pattern predicted by perturbative QCD had 
not been compared systematically with data. At PEP energies the TPC col- 
laboration [Hj found that their data - though in qualitative agreement with 
QCD expectations - showed some quantitative differences with the asymp- 
totic analytical results, which could have possibly been caused by the low 
purity of the quark/gluon identification of the jets. At LEP there has been 
considerable progress in jet identification, in particular by DELPHI [9]; their 
result for the "Mercedes configuration" has been compared successfully [3] 
with the analytic calculation Comparison with the limited results by 

OPAL averaged over a large range of jet angles, has been moderately 
successful. 

In this paper we compare the experimental results for the various sym- 
metric and asymmetric angular configurations of qqg presented by DELPHI 
and OPAL for the full angular region, along with the results for qq'j and two 
jet events with theoretical predictions. Furthermore, we discuss the energy 
and jet resolution dependence of the effect which has been ignored in the 
previous analyses. 
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2 Angular pattern of particle flows 



2.1 Perturbative results and dependence on the jet se- 
lection 

Let us begin by recalling the main ideas in the perturbatively-based cal- 
culation of particle flows in e + e~ — > qqg at high energies in comparison with 
e + e~ — > qqj. We consider the configurations where all the angles 0jj be- 
tween jets are large (i — {-I — 1} = {ggg} or i = {H — } = {qq})- Within 
the perturbative picture, the angular distribution of soft inter-jet hadrons 
is calculated from the distribution of soft gluons radiated coherently off the 
colour antenna formed by the primary emitters (q,q and g) or (q and q). 

The angular distribution of a secondary soft gluon, g 2 , is derived in low- 
est order perturbation theory from the corresponding Feynman diagrams of 
single gluon emission off the primary partons. For the qq^ process one finds 
neglecting the recoil 
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with 4- momenta p+,P- and k 2 and soft gluon energy and transverse momen- 
tum E 2 and k t . 

A modification is necessary to take into account the fact that the "de- 
tected" gluon belongs to a parton jet. The results of the calculation for qqj 
can be written as 
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Here the angular distribution of soft gluon bremsstrahlung from the "antenna" 
dipole (ij) is obtained from and reads 
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This angular distribution is the same as that for photon bremsstrahlung in 
QED - as it occurs, for example, in e + e~ pair creation. Unlike the QED case, 
however, there are also the QCD-specific colour factors and the "cascading 
factor" N' g (Y), which takes into account the fact that g 2 is part of a jet. It 
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represents the derivative of particle multiplicity in a gluon jet with respect to 
Y = Iq^Kt/Qo) at the appropriate maximal transverse momentum scale 
and cut-off Qq. The parameter Qo is determined from a fit to the multiplicity 
data together with the QCD scale A (Q > A, see, for example 

The scale Kt depends on the way the jets are selected. To see this, we 
write down the perturbative expansion for the emission of gluon g 2 in (J2J) 
from the primary parton p as follows (see also [T2]) 

(5) 

The first term corresponds to the Born result for direct gluon emission in 
the second one to the emission through the intermediate parton a. We 
denote by and Q a b the transverse momentum and angle respectively of 
parton a with respect to parton b, in addition 6 ap = 9 a ; the leading order 
multiplicity anomalous dimension is given by 7o(&t) = ^N C 0i s {k t ) /it. Eq. ((SJ) 
is written in DLA, but it can also be generalized to MLLA after inserting 
the appropriate splitting functions in all intermediate branchings. 

In the second equation in (J5J) we have replaced the integral over the angle 
@a P by the integral over Q a2 , as the leading contribution to the a integral 
comes from the region of quasi-collinear emission n a ~ n 2 ; then we also 
replace the corresponding angles involving particle a in Furthermore, 
the order of integration between E a and E 2 is interchanged (E 2 < E a ). Now, 
the angular integral in the second term has to satisfy the "angular ordering" 
requirement, i.e. the emission angle of the intermediate parton a should 
respect 

@2a < @a P ; therefore, also Q ap > <d 2p /2. (6) 

Then the angular integral can be performed with Q a2 < Q 2 in pole approxi- 
mation. 1 

Next we discuss the limits of the energy integrations in (JHJ)- The lower 
limit E min is determined by the transverse momentum cut-off. For gluon 
emission inside a well separated jet one requires k t > Qo, therefore E min ps 
Qo I @2 in the small angle approximation. In case of emission from a boosted 

lr The discussion of boundaries in these integrals is the same as in the derivation of the 
inclusive spectra, see for example, Ref. |LS| . 
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qq "dipole" as in (JTJ) an azimuthal dependence of the effective cut-off Qq 
around the jet axis is generated. To estimate this effect we consider a gener- 
alization of the cut-off restriction following from Q 

p _ 2(p + k)(p.k) _ 2E 2 2 (f-cose +2 )(l-cos9_ 2 ) 2 
' ~ ~~P^V- ~ f - cos + _ > Qo ' {V 

which approaches the limit k t = k t for the qq back-to-back configuration. 
For our applications the effective k t cut-off Qo will then change by up to 
a factor 2. In the asymptotic expansion of multiplicity |14f2j In N g (Y) = 
c\\/Y + C2 In Y + C3/ \/Y such a rescaling Q — > /3Qo would modify only the 
C3 term, therefore the leading (DLA) and next-to-leading (MLLA) results 
remain unchanged. Numerically, for the present analysis we estimate the 
effect on N' in p|). ([5j) from the variation of Qo to be at the 10-20% level 
where we take into account the data in Fig. 4a below. In our application 
of analytic high energy approximations (large Kt/Qo) we continue therefore 
with the cut-off k t > Q in the following. 

The upper limit E max of the energy integral depends on the jet selection 
procedure and we consider two cases 
a) no momentum restriction in jet selection 

In the simplest case, we require 2 or 3 respectively energetic particles in 
the given angular directions which define the jet directions. A more precise 
definition is possible using the energy-energy-multiplicity or energy-energy- 
energy- multiplicity correlations for 2 and 3 jets |15|3j . For two jets one 
considers all pairs of particles (ij) weighted with their energies, and for a given 
relative angle Oy, one studies the angular distribution of all soft particles 
in the event. In the same way one proceeds for 3-jet events with triples of 
energy weighted particles and two relative angles. 

In these cases there is no specific restriction on the energy of the triggered 
soft particles whose angular distribution is investigated. The maximum en- 
ergy E max of these particles is then of order of the respective jet energies 
Ejet, say y/s/2 or \/s/3 in 2 and 3 jet events at total cms energy y/s of the 
jet system, and the same applies to the intermediate partons a. 

b) selection of jets for a given resolution parameter y cut 
We consider here the "Durham algorithm" ^Hj which selects the jets accord- 
ing to a predefined minimal relative transverse momentum, approximately 
Kt > \fUcutS- In this case, the y cu t selection restricts the transverse momenta 
of emitted particles as well as the transverse momenta of the intermedi- 
ate partons a in © in the same way (fc* p < K T ), whereas the longitudinal 
momenta of the latter partons a are only limited by 0{\fs). However, an 
additional restriction comes from the angular ordering (j§J, which is violated 
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for intermediate partons a with large energies and, thus, small angles <d ap . In 
our logarithmic approximation we take the same upper limit for the E a and 
E 2 integrals. Similarly, the angular integral is limited by a 2 < ©2 within the 
same accuracy. Then, the second equation in J5J) can be resummed, resulting 
in the multiplicity N 9 (Kt /Qq) 

87rdN a a-y 1 , ^. f Emax dE 2 9/W Wr /W ,~ X 

_^ „ _ H / + _(rr 2 ) Jf^ -^7„ 2 (^)iV 9 (*yQ„). (8) 

This integral can be expressed in terms of N'(Y) using the evolution equation 
for multiplicity, which then leads to Eq. Q. We note that this replacement 
is also possible in MLLA. 2 

The upper limit of integration E max is either given by the maximal energy 
Ej et of the jet or depends on the jet resolution parameter y cut 

case a (no restriction) Y = ha.(Ej et Q m /Q ) (9) 

case b (y cut restriction) Y = H E ^/Qo) © m < ^E/E jet 

\l^{y/y C uts/Qo) otherwise 

Here 9 m is the angle between the soft gluon g 2 and the jet closest in angle. 
As the multiplicity N g depends on the maximum transverse momentum scale, 
the Durham Kt - algorithm naturally provides simple results. Note that in 
case of y cut selection (case b), the N'(Y) factor becomes independent on the 
emission angle O m away from the jet directions, and, therefore, the soft gluon 
angular distribution is given simply by the Born-term factors W{n2) in (|2|Hjl . 

2.2 Projection onto the event plane 

We are now interested in the projection of the soft gluon g 2 angular distri- 
bution onto the event plane defined by the momentum vectors of the qqg{^) 
systems. The azimuthal angle <fi in the plane is defined to be zero in the 
direction of q, becoming positive in the direction of q for definiteness (in 
our applications q and q are not distinguishable), see Fig. [TJ The projection 
of the particle density onto the event plane is obtained in the case of qq'j 
from the integral over cos 02, where 02 is the polar angle with respect to the 

2 The MLLA evolution equation [1 1 1 is usually written in small angle approximation 
with Kt = K@, the continuation to larger angles is not unambiguous. 
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Figure 1: Definition of azimuthal angles in the event plane of three jet event. 
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the angles a, j3 are given below in terms of 4>. 

Replacing 7 by the hard gluon gi one obtains an additional particle flow 
in the g\ direction, and from the integration of Eq. Q over cos02 one derives 



W- 



±11 



N, 



c 



a+iV(a, 7) + a-iV(P, 7) - a ~\ 



N, 



c 



(12) 



For an arbitrary choice of the angles a,/3, 7 in the projection formulae 
()ll |) . (jl2J) have to be defined as minimal angles between the directions of the 
soft gluon (#2) and the directions of q, q and g respectively, with all angles 
within the interval [0, tt]. We find the following definition of the angles a, (3, 7 
to satisfy this requirement and to apply in all angular sectors with cb running 
from to 2tt (see also Fig. [TJ) 
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Eq. (|12j) corresponds to an incoherent superposition of two Lorentz- 
boosted dipoles, one between the gluon and the quark and one between the 
gluon and the anti-quark, modified by the negative correction of 0(1/ N^). 
This expression shows a depletion opposite to the gluon direction. 
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2.3 Comparison with experiment 

Let us now compare these analytical formulae with experimental data. The 
DELPHI collaboration [Oj has presented results on three jet events with two 
identified 6-quark jets allowing for a high gluon jet purity (~ 94%). Sym- 
metric angular configurations have been selected - either "Mercedes" or "Y- 
symmetric" events with angles confined to an angular region of ±15°. Fur- 
thermore, the data on the final state qcpf in the Y-configuration have been 
presented. These well defined event classes are compared in Fig. El with the 
above analytical formulae. 

The angle <fi = is in the direction of the most energetic quark jet. For 
Mercedes events, we chose the other jet angles at the observed peak positions 
- slightly different (< 8°) from the symmetric values. The events had been 
selected with fixed y cu t — 0.01 according to the Durham algorithm. As a 
result, we can neglect the variation of the cascading factor N' in (J2J) and (JHJ 
(for angles > 20° away from the jet directions). The shape of the angular 
distribution is then predicted without a free parameter; in Fig. El only one 
overall normalization factor has been fitted in the description of the three 
distributions. 

The general characteristics of the distributions are well reproduced, in 
particular the different heights of the inter-jet valleys. In addition, the larger 
width of the gluon jet in comparison with the central quark jet width becomes 
visible. The quark-jets in a) and b) are a bit broader than perturbatively 
predicted, as can be expected from the b-quark jet selection; the quark-jets 
at the larger angle <fi in a) and b) are even broader because of the angular 
fluctuation in the selection of the jets. This successfull description nicely 
demonstrates the similarity of particle flow in "strong interaction" processes 
to photon flow in electromagnetic interactions - both are caused by the un- 
derlying gauge boson bremsstrahlung in QED and QCD. In c) the formulae 
are equivalent, in a) and b) they are modified by the different colour charges 
for quarks and gluons - a degree of freedom not available in QED. An interest- 
ing QCD interference effect is the negative 1/Nq term in Unfortunately, 
according to our analysis, the presented data are not sufficient for a reliable 
observation of this term. 

A comparison of qqg and qcpj distributions between the quark and anti- 
quark jets has been presented by OPAL 10 . In the analysis of qqg events, 
the jets are ordered according to their energies. The jet with the lowest 
energy in the qqg system is taken as the gluon jet and two intervals of this 
energy E 3 are selected around E 3 = 10 and E 3 = 20 GeV, furthermore, 
© + _ ~ 165° from which we can estimate 0i_ and Oi+. In the calculation we 
take into account the estimated purities of gluon jets, p g , by superimposing 
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b) qqg events with Y topology 
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Figure 2: Charged particle flow within various multi-jet configurations as 
measured by DELPHI [9 j (a) "Mercedes" qqg events, (b) "Y-symmetric" qqg 
events and c) qqj events, in comparison with analytical QCD predictions; the 
curves correspond to the lowest order QCD soft bremsstrahlung formulae. 
The relative angles between the jets are taken in a) as O i+ = 125°, Q + _ = 
122.5° and 0i_ = 112.5°; in b) and c) as 0+_ = 1+ = 150°. 



the distribution (J2J) with the one obtained by swapping the second and third 

jet, i.e. by swapping H and 1+ and changing <fi to 2ir — <ft. The original 

distribution gets the weight p g = 0.92 for the sample with E 3 = 10 and 
p g = 0.74 for E 3 = 20 GeV, according to the reported purities. 

In Fig. El we show the distribution in the rescaled azimuthal angle X = 

0/0^ for the two processes as well as for their ratios in the two intervals of 

E 3 . The different heights of both distributions near X ~ 0.5 is rather well 
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reproduced, as is the variation with energy E 3 and jet angle 0i_. In partic- 
ular, there is a sizeable asymmetry in the distribution which would be even 
more pronounced if the gluon jets were identified with higher purity. Since 
the data in both intervals still average a sizable range of different angular 
configurations, some deviations from the predictions are to be expected. We 
therefore allowed an increase of relative normalization by 15% of the curves 
in Fig El b) over Fig |3] a) to improve the description. The predictions are 
seen to fall below the data for dn/dX near the jet directions, which may be 
due to the averaging over the jets with different energies and angular widths. 
On the other hand, the difference in the splitting of the two distributions, in 
their asymmetry, and in their absolute height for the two angular configura- 
tions in a) and b) are well reproduced. The effects related to configuration 
asymmetry have been studied here for the first time; these effects could be 
investigated more precisely in the future if smaller intervals in the relative 
jet angles were chosen. 

3 Energy dependence of particle angular flow 

The energy (\fs) and y cut dependences are determined by the cascading factor 
N'(Y) according to ljnjl. (|IU}) . Angular distributions for particles produced in 
the hemispheres of quark and gluon jets without y cut restriction are available 
at different energies from e + e~ — > qq and at LEP-1 energy from e + e~ — > bbg. 
Instead of the distribution in angle we may also study the distribution in 
pseudo-rapidity y = — lntan(0/2). For quark and gluon jets (A = q,g) we 
have 



This relation is not limited to the leading order result N q /N g = Cf/Ca as in 
Eq. (J2J) but allows for higher order logarithmic effects. Relation (fT3|) follows 
from the MLLA result that the multiplicity depends only on the maximum 
Kt in the jet, i.e. on the product Ej et Q in the small angle approximation. 3 
Then, the variation of multiplicity with y ~ — ln(0/2) is the same as the 
one with lnEj et ~ y ma x = ln(2£y et /m) or, equivalently, the particle density 
depends on rapidity and energy only through the scaling variable 



3 Note that the remarkable MLLA prediction |2I."{| of such Ej et Q scaling behaviour 
has been recently well confirmed experimentally by the CDF collaboration JJj in the 
measurements of inclusive charged particle production in restricted cones around the jet 
direction. 
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Figure 3: Charged particle flow in between the two quark jets of the qqg and 
qqj final states, as measured by OPAL jTUj, for two samples of energies E 3 of 
the lowest momentum jet together with the corresponding ratios. The curves 
represent the lowest order QCD bremsstrahlung formulae for two intervals of 
the lowest energy E 3 of the jets; the overall normalization has been adjusted; 
in addition, the curves in b) are increased by 15% (see text). In our calcula- 
tions we use for the first interval E% = 10 GeV, + _ = 165°, 0]_ = 67° and 
6i+ = 128°, and for the second one E 3 = 20 GeV, 6+_ = 165°, 6i_ = 35° 
and 9i+ = 160°. 



where we set \i = 1 GeV and yfs = 2Ej et . Eq. |T5|) predicts both the energy 
dependence at fixed angle (rapidity) and the angular dependence at fixed 
energy. Recall that within the MLLA approach |2fllj the multiplicity N de- 
pends only on the phenomenological parameter Q ~ A and a normalization 
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Figure 4: Distribution in rapidity y (a) at fixed rapidity y = 2 versus 
energy and (b) as a function of the shifted rapidity x = y — In (^/s/ GeV) 
at different energies y 7 ? (assuming pion mass for all charged particles), in 
comparison with solutions of the MLLA evolution equations for N' by Dremin 
and Gary (DG) [H] and by Lupia and Ochs (LO) [EH]; data by TASSO j20j 
and OPAL EQ. 

factor, then the angular distribution is uniquely predicted by (|T5|). 

These predictions have been tested by comparing the QCD results for 
N' A with data from TASSO and OPAL (here from w<is-quark jets only to 
avoid the heavy-quark contribution at the Z mass) in Fig. 0J The curves 
in Fig. |3] represent the high energy QCD 3NLO prediction (DG JSj) 4 and 
the numerical solution (LO ^H]) of the full MLLA evolution equations for 
quark and gluon jets ^1] . In both calculations the parameters have been 

4 we selected the fit with parameters nj—3, A=0.67 GeV and K=0.288 for a gluon jet 
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determined from a fit to e + e~ data, and the two results agree for quark jets 
within ~ 3% in the energy region considered (10 < ^/s < 200 GeV). The 
data in Fig. show the predicted increase with energy. In Fig. HJd we show 
the comparison for the rapidity distributions at different energies using the 
scaling variable x in (jlfij) . At each energy, the lower bound of the rapidity 
interval is y — 0.8. Above this value the scaling property is approximately 
satisfied, and, again, a good description by the theoretical calculation of 
N'(y) is obtained in absolute terms. However, some discrepency at low y 
is seen for the OPAL data. 5 Approximate agreement is also obtained with 
the distribution in gluon jets from OPAL [21] in the given rapidity range. 
For larger rapidities, the distibution would fall below the distribution for 
quarks as one may expect from energy conservation and recoil effects given 
the increased central production of particles. The curve shown corresponds 
to the common solution of the evolution equations for quark and gluon jets 
(LO PEj); the solution DG [IB] with the quark jet parameters would result 
in a prediction larger by 20%; with readjusted parameters a good description 
is obtained again. 

We conclude that the energy dependence of the angular distribution, given 
by N'(Y) in absolute terms, is reasonably well reproduced. A further im- 
provement would require the treatment of (i) large angle and recoil correc- 
tions, (ii) dependence on jet axis definition (iii) heavy quark contribution 
and hadron mass effects. 

The energy dependent effects discussed here are also expected in multi- 
jet events as in Eqs. (JSJEJ). Different analysis methods and configurations 
selected by different groups do not allow a meaningful comparison yet. An- 
other interesting test will be the dependence of angular distributions on the 
jet resolution parameter y cut in (JTUj). 

4 Conclusion 

The angular distribution of particles in two and three jet events follows closely 
the expectations based on the lowest order bremsstrahlung formulae of QCD. 
These formulae can be represented as sums over dipole radiation terms (@J), 
the same as in QED, modified by the proper colour factors for quarks and 
gluons. In addition, there is an effect of parton branching which modifies 
these formulae by the "cascading factor" N'(Y), the logarithmic derivative 
of the multiplicity. For fixed resolution parameter y cut , this factor does not 

5 OPAL defines rapidity with respect to the sphericity axis, while TASSO defines ra- 
pidity with respect to the thrust axis; without the uds selection and using the thrust axis 
to define rapidity ALEPH 22 data would lie above the curves. 
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modify the basic angular distribution in a wide range. On the other hand, 
it implies an additional dependence on energy and/or jet resolution. The 
predicted dependence on N' has been clearly established for two-jet events. 
These results demonstrate that soft multi-hadron production follows the sim- 
ple expectations of perturbative QCD bremsstrahlung and confirm the close 
connection between the colour flows and observable flows of hadrons. 
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